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Abstract: We study in detail the scalar-field Hawking radiation emitted into the bulk by a 
higher-dimensional, rotating black hole. We numerically compute the angular eigenvalues, 
and solve the radial equation of motion in order to find transmission factors. The latter 
are found to be enhanced by the angular momentum of the black hole, and to exhibit the 
well-known effect of superradiance. The corresponding power spectra for scalar fields show 
an enhancement with the number of dimensions, as in the non-rotating case. On the other 
hand, the proportion of the total (i.e., bulk-l-brane) power that is emitted into the bulk 
decreases monotonically with the angular momentum. We compute the total mass loss 
rate of the black hole for a variety of black- hole angular momenta and bulk dimensions, 
and find that, in all cases, the bulk emission remains significantly smaller than the brane 
emission. The angular-momentum loss rate is also computed and found to have a smaller 
value in the bulk than on the brane. 



Keywords: Large Extra Dimensions, Black Holes, Beyond Standard Model. 



Contents 



1. Introduction 1 

2. Theoretical background 3 

2.1 Scalar field equations 4 

2.2 Hawking emission 5 

3. Numerical methods 9 

3.1 Angular eigenvalues 9 

3.2 Transmission factors 10 

4. Results 11 

4.1 Non-rotating black holes (a^ = 0) 11 

4.2 Rotating black holes, > 12 

4.2.1 Transmission factors 13 

4.2.2 Power spectra 15 

4.2.3 Angular momentum spectra 19 

4.3 Consistency with other studies 22 

5. Discussion and conclusions 23 



1. Introduction 

Theories with Large Extra Dimensions [1, 2] have proven attractive to theorists in recent 
years. Of particular interest are models in which matter is confined to a four-dimensional 
hypersurface (the brane), but gravity (and possibly scalar fields) is free to propagate in 
a higher-dimensional compact space (the bulk). Phenomenologically, such models are in- 
triguing because they appear to resolve the so-called hierarchy problem. In other words, 
they explain why gravity is observed to be so much weaker (on the macroscale) than the 
other forces. 

In these models, the fundamental energy scale of gravity {M^,) is related to the Planck 
energy {Mp) by M^+" ~ MpR~'"', where R and n are the size and number of extra 
dimensions. Hence, M* may be many orders of magnitude lower than the Planck energy 
Mp. This raises the tantalising prospect that black holes may be created through 'trans- 
Planckian' particle collisions [3]. It has been suggested that may be probed by high- 
energy cosmic ray collisions [4], or the next generation of particle accelerators [5]. If 
created, mini black holes would evaporate rapidly through emitting Hawking radiation 
[6]. Experimental detection of Hawking emission would provide a clear signal of black 
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hole creation. Aeeurate measurement of the power spectrum would enable us to deduce 
the properties of the underlying spacetime itself. For these reasons, accurate theoretical 
modelling of the emission spectrum has become a high priority, and it has received much 
attention in recent years. 

Black holes are thought to undergo four stages of evaporation: the so-called 'balding', 
'spin-down', 'Schwarzschild', and 'Planck' phases [7, 8]. Thanks to a number of analytical 
[9, 10, 11] and numerical [12, 13] studies, emission from the Schwarzschild phase is now well 
described. Attention has recently shifted to emission in the rotating phase [14, 15, 16, 17]. 

In this paper, we conduct a numerical study of scalar emission in the bulk by a ro- 
tating higher-dimensional black hole. Our work complements a range of existing studies 
of the 'spin-down' phase of black hole evolution. In previous papers in this series we have 
examined on-thc-brane emission of scalars [15], fcrmions [16] and photons [17] by a rotating 
black hole. In parallel, other groups have pursued complementary lines of inquiry [18]. The 
bulk emission of scalars from the rotating phase has received increasing attention. Analytic 
approximations for the greybody factors (transmission factors) were recently derived [19]. 
The analytic work is complemented by two recent numerical studies [20, 21] of bulk scalar 
emission from 5D and 6D rotating holes. Here, we present comprehensive new exact results 
for bulk dimensionalities n = 1, 2, ... 6, and a range of black hole angular momenta. 

This study aims to investigate the claim that black holes radiate mainly on the brane 
[22]. Studies of the Schwarzschild phase of evolution found that graviton [13] and scalar 
[12] emission in the bulk is small compared with emission of standard model fields on the 
brane, due chiefly to the multiplicity of standard model particle species. Hence, only a 
small fraction of the black hole energy is 'lost' in the bulk. It is widely suspected, but 
unproven, that the same remains true when the black hole is rotating. As the amount of 
energy emitted in the bulk inevitably defines the one on the brane - and its observable 
signatures, this point should be duly clarified. Two species of particles are usually assumed 
to propagate in the bulk, gravitons and scalars. Whilst it is possible [23] to derive a set of 
partial differential equations (PDEs) describing the gravitational perturbations of a higher- 
dimensional rotating black hole, it is not yet clear whether the PDEs can be decoupled 
into ordinary differential equations (ODEs) for all modes (specifically, the scalar and vector 
modes) when the black hole has a single angular momentum parameter^. On the other 
hand, in the case of scalars, the decoupling of PDEs is indeed realised, and a comprehensive 
study of their emission can be performed. Apart from the importance in its own right of 
the scalar emission, the study of the scalar field in the bulk can also offer a qualitative 
understanding of higher-dimensional graviton emission from a rotating black hole. In the 
sections that follow, we compute the strength of scalar emission in brane and bulk channels 
in the rotating phase, and calculate the proportion of the scalar energy which enters the 
bulk. We make the assumptions of a minimally-coupled field, negligible brane tension, and 
a black hole with just one non-zero angular momentum component. 

The paper is structured as follows. In section 2 we outline the relevant theory, viz. a 
mathematical description of a rotating black hole in a bulk spacetime; the decomposition 

^In the case of equal angular momenta, the ODEs have been separated for some modes - see [24] and 
[25]. 
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of the scalar field into radial and angular parts; and the relevant formulae for Hawking 
emission. In section 3 the numerical methods employed to calculate angular eigenvalues 
and transmission factors are described. In section 4 we present our numerical results. We 
start with a brief review of the Schwarzschild phase results, before moving on to consider 
the rotating phase. The effect of superradiance on the transmission factors is examined. 
We present branc and bulk power spectra for a range of n and black hole angular momenta 
a. We compute the total scalar emission rate, and the fraction which is radiated into the 
bulk. We also consider the loss of angular momentum from the spinning hole. Finally, in 
section 5 we conclude with a discussion of the physical implications of our work. 

2. Theoretical background 

The well-known Myers-Perry solution [26] describes the gravitational field of a (4 -|- n)- 
dimensional uncharged rotating black hole. We will confine our attention to black holes 
created by the collision of particles on the brane. Hence we assume that the black hole 
metric has only one non-zero angular momentum component, in a plane parallel to the 
brane. The line-element takes the form 

o / u, \ r. 2an sin^ 9 So q 

ds'=-{l- ^) dt' - -^-^dt # + -dr' + S de' 

+ (^r^ + 0" + ^^^^) sin2 9 d^^ + cos^ 6 d^l (2.1) 

where 

A = r2 + a2--^, S = r2 + a2cos20. (2.2) 

Here, 

dVLl{Q^,92,...,9r,-^A) = 

d9l_^ + sin^ On-i (del_2 + sin2 9n-2 ( • • • + sin^ 62 ( dOj + sin 9ld^'^) . . . ) ^ (2.3) 

is the line-element on a unit n-sphere. In equation (2.1), the angle ijj is the azimuthal 
angle round the axis of rotation of the black hole in the brane (for four-dimensional black 
holes, this angle is usually denoted cf)). The mass Mbh of the black hole and its angular 
momentum J are proportional to /x and afj,, respectively, 

Mbh = m^/x, J = -^Mnna, (2.4) 
IottG n + 2 

where An+2 = 27r("+3)/Vr[(n + 3)/2] is the area of an (n -|- 2)-dimensional unit sphere, 
and G is the (4 -|- n) -dimensional version of Newton's constant. Equation (2.4) effectively 
defines the parameters /i and a which appear in the metric (2.1). 

The radius of the black hole's event horizon, r/j, is the largest, positive root of A(r) = 0. 
The horizon radius can be employed to define a useful set of dimensionless variables, 
which we denote with the subscript d, 

rd = r/rh, = a/rh, ujd = ^rh- (2.5) 
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For n > 1, there is one unique root in the region r > 0, and it can be obtained from the 
implicit equation 

For n > 1 this root exists for any value of Ud- However, it has been shown [27] that the 
maximum value of for black holes created by particle collision is = n/2 + 1. 

2.1 Scalar field equations 

In this paper, we study the Hawking emission of scalar fields into the higher-dimensional 
'bulk' spacetime. Therefore we must first consider the equation of motion for a massless 
scalar field $ propagating in the bulk, with minimal coupling to the geometry, which 
satisfies the field equation: 

^ d^{V^g^''d,^)=0. (2.7) 



V-9 

Here, ^/—g is the determinant of the metric. 



n-l 



^= Sr"sin^cos"^ JJsin^^j. (2.8) 



i=l 



Substituting in (2.7) the factorization ansatz 

$^ = e-"^'e'"'^RA{r)SAi9)Yjni9i, . . . , 9n-u 4>), (2-9) 
where A = {Z, m,j, w}, yields three coupled second-order equations. The first equation. 



n— 1 ^ 

gn-^sin^e.^^^ 



o^<t>yj 



sm 



^ rr^-T / +.^-(i+n-l)y.n = 0, (2.10) 



defines the hyperspherical harmonics on the n-sphere, Yjn [28]. The remaining two equa- 
tions determine the angular Sa{9) and radial RA{r) mode functions. They are [19]: 

-—L^Oeisin 9 cos" 9 5,5a) + (io^a'' cos^ 9 - ^ - ^^^^"l'^^ + Ea] Sa = 0, 
sm 6 cos" 9 \ sm 9 cos^ 9 J 

(2.11) 



—dr (r"A5^i?A) + (^- ■^'(■^' + " ^)"^ _ \Ea + 2amuj - a^w^l ^ i?A = 0, (2.12) 



and 

^,-dr{r-/^drRA)+{-^ 

where K = {r'^ + a^)a; — am. The angular (2.11) and radial (2.12) equations are coupled 
through the angular eigenvalue Ea, which depends on the spheroidicity parameter auj. 
There is no known closed analytic form for the eigenvalues, but at least three methods 
are available for their calculation. First, and simplest, Ea may be expressed as a power 
scries in auj. The coefficients of the power series up to 4th order in auj are given in [29]. 
Alternative numerical methods are Leaver's method of continued fractions [30], and a 
numerical matching method, described fully in section 3.1. 
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Note that the angular equation (2.11) together with the boundary conditions of reg- 
ularity of the solution at the points ^ = and ^ = 7r/2 constitutes a Sturm-Liouville 
eigenvalue problem. Orthogonality of the solutions with the appropriate measure follows, 
and we normalize them so that 

/ desmecOs'^eSyrnjojSlmju, = Sl'l. (2.13) 

Jo 

Similarly, the hyper spherical harmonics are orthogonal and we normalize them as 

/n— 1 n— 1 
# n n si""" ^kYfnY*^ = 5yj. (2.14) 
i=l fe=l 

2.2 Hawking emission 

The quantization of the scalar field in the higher-dimensional curved space-time follows 
through in essentially the same manner as in the 4-D case. We will therefore merely 
indicate the main steps in the quantization and derivation of the flux formulae and we 
refer the interested reader to, for example, papers [31, 32, 33] for the 4-D case and to [14] 
for the n = 1 case. 

We define the inner product of two solutions ^ and $ of the scalar equation (2.7) in 
the bulk as 

^) = \j^ d^^V^ i^*^^ - ^*^^) , (2.15) 

where a semicolon denotes the covariant derivative, S is a complete Cauchy surface and 
dS'* is a future-directed normal to E (in this section only, the symbol S will be used to 
denote a complete Cauchy surface rather than the function in equation (2.2)). Using the 
scalar field equation it can be checked that this inner product is independent of the choice 
of the complete Cauchy surface S. 

We will follow Myers and Perry [26], and define null coordinates u and v as 

u = t — r*, ■u = t-|-r*, (2-16) 

where r* is a "tortoise" radial co-ordinate defined as 



and K is the surface gravity, which is given by [26], 

^^(n+l) + (n-l)a^^^.^r A 



2(1 -I- a2)rh r-^^h 



2{rl + a?){r-rh). 



(2.18) 



The null congruences are tangent to surfaces of constant u or i; at the horizon {r = Vh) 
and at infinity (r — oo), but this is generally not the case for other values of r. We can 
also define the Kruskal-like co-ordinates, 

U = -e""'^, V = e""^. (2.19) 
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We will choose E = U as the Cauchy surface, where is the past null infinity 
given by ?7 = — oo and is the past horizon given by y = 0. The surface element in 
these regions can be calculated, 



where I and n are two null vector fields, 

d (r^ + c? 



^ "o^^^ o It^ ^^^^ # dQx--- dOn-i d4> at 7^- 
(r^ + cos^ 0)(7 

^ nf'dVd9di;d9i...den-i( 



A 

(r2 + a2) 



at J" 



(9 a did 

dt (r^ + a^) (J^/jJ 9r ' 







d_ 



5^ 
dr 



(2.20a) 
(2.20b) 

(2.21a) 
(2.21b) 



The radial equation (2.12) has singular points at r = r^j and infinity. We can define 
two sets of linearly-independent radial solutions (valid for n > 0) with the asymptotic 
behaviours, 



■(rd-l)-^'^/2«, 



1 

l+n/2 



+00, 



''A 



B'^ivd - 1)-'^/^'^ + B'^\rd - l)+^'^/2'^, 



1 „+iw^r^ 
l+n/2 ^ ) 



+c>o, 



(2.22a) 



(2.22b) 



where AJ^, -B](^ and B^*^* are complex constants, Hi = oo — mQh and Cl^ = a/{r^ + a^) 

is the angular speed of the horizon. We define the 'in' solutions for w > and the 'up' 
solutions for cD > 0. 

The field mode solutions $™/"p constructed via equation (2.9) from equations (2.22) 
satisfy the normalization conditions 

= 0, 



(2.23) 



after an appropriate normalization factor is included, where symbols • and •' denote either 
the 'in' or the 'up' solutions. The normalization factors can be calculated using the surface 
elements (2.20) and, when included, the field modes are explicitly given by 
1 



^7 



27rV^r"/'+V 



e-'-'e'^^ir;^{r)SA{e)Yjn{ei, en-i,<p), 



(2.24) 
(2.25) 



where the orthogonality relations (2.13) and (2.14) have been used. 

The field mode solutions 'in' and 'up' form a complete set of orthonornial solutions to 
the scalar equation (2.7) in the outer region of the space-time. A general solution of the 
scalar equation can therefore be written as a linear combination of this basis 



A •=in,up 



(2.26) 
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where the integration over the frequency is for uj > for the 'in' modes and for uj > for 
the 'up' modes. Promoting the fields and the Fourier coefficients a\ to operators a\, using 
the normahzation conditions (2.23) and imposing the standard canonical commutation 
relations for the field and its canonical momenta leads to the commutation relations 



0, 



(2.27) 



In the eternal version of the space-time (i.e. in the full analytic extension of the Myers- 
Perry metric) we can construct from the field modes ^^^^"^ a complete orthonormal basis 
of modes which are positive- frequency with respect to both the affine parameter U on the 
past horizon 7i~ and the affine parameter t on past null infinity J~ . The past Unruh 
state \U~) is defined [31, 34] as the quantum state which is annihilated by the Fourier 
coefficients - as operators - of this basis of modes. This state models an evaporating black 
hole and is therefore the relevant one for computing the Hawking emission. 
The stress-energy tensor for a minimally coupled scalar field is given by 



1 



(2.28) 



It can be proven [31, 32] that the expectation value of the stress-energy tensor operator 
when the scalar field is in the past Unruh state is given by 



u- 



u- 




duj coth 



2Th 



' ^A 



+ 



f 

JO 



duT, 



IJ,V 



with the Hawking temperature 



(2.29) 



(2.30) 



It is well-known that the expectation value of the stress-energy tensor is generally divergent 
and it must be renormalized appropriately. One method of renormalization is the technique 
of covariant geodesic point separation. In this approach, the pair of fields appearing in each 
quadratic term in the stress-energy tensor (2.28) are evaluated at separate points, x and 
x' . This defines a bi-tensor T^y{x,x'). The expectation value of this bi-tensor is found for 
the desired quantum state, and certain purely geometric terms T^"{x,x') are subtracted 
from it. Finally, the two points are moved together again, i.e. the coincidence limit 
x' ^ X is taken. In 4-D, Christensen [35, 36] calculated the geometric subtraction terms 
using the Schwinger-DeWitt expansion [37, 38] for the Hadamard elementary function. 
The Schwinger-DeWitt expansion for the non-rotating, higher-dimensional case is given in 
[39]. In principle this could be used to compute the geometric subtraction terms in higher 
dimensions, following [35, 36], but fortunately we do not require their precise form in our 
analysis. 

In [33] it was shown that those components of the unrenormalized expectation vahic 
of the stress-energy tensor with one index either t or %1j and the other index r (as in (2.31) 
below) are not divergent and so do not require renormalization. The proof is based on two 
properties. Firstly, we use the symmetry of the metric under (t, ■i/') (— —i^)- Secondly, 
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each of the geometric subtraction terms contains an even number of covariant derivatives cj^ 
of the biscalar of geodetic interval [40] after an average is taken over a separation in the 
and — cr'^ directions. By choosing a point-sphtting in the radial direction, these two proper- 
ties then ensure that the geometric subtraction terms T^^^{x, x') and T^^{x, x') vanish (see 
[33] for the full argument). Even without explicitly computing the geometric subtraction 
terms, it can be shown that these two properties remain true in higher dimensions, so that 
the argument of [33] follows through, and we have T^"{x,x') = T^^{x,x') = independent 
of n, the number of extra dimensions. 

Having established that the tr and ij^r components do not require renormalization, the 
fluxes of energy and angular momentum are given by 

'^ = A jxi^^2{u-\f;f\u-). (2.31) 

where S is any surface of constant t and r. The total rate of energy emission of massless 
scalars into the bulk is found by taking a sum over all angular modes. Taking for the 
surface S in (2.31) a (2 + n)-sphere boundary at infinity and using equations (2.29) and 
(2.22), it can be shown [14, 19] that 

<fE ^^_y (2.32) 



'-tr 



u- 



dtduj 2n ^ expfcj/Tij) — 1 

and the emission rate of angular momentum into the bulk is given by 

^ J 1 X ^ m 

dtd^^2^^ exp{Co/TH) - 1 ' ^ ' 

l,j,m 

Here, 7a is the energy-dependent transmission factor (also known as 'greybody factor') 
and is given by the ratio of the ingoing and outgoing fluxes, 

2 



rA = i 



,iout 



4 in 



(2.34) 



In addition, A'^ is a degeneracy factor (given in equation (11) of [28]) which accounts for 
the multiplicity of modes in the bulk, 

^ (2j + n-l)(i + «-2). 

' jHn-iy. ' ' 

All parts of equation (2.32) can be determined analytically, except for the transmission 
factors 7a • Qualitatively, the transmission factor is a dimensionless measure of the propor- 
tion of a given mode that can escape from the black hole horizon to reach inflnity. Hence, 
7a = 1 corresponds to 'total' emission, and 7a = to 'total' back-reflection of Hawking 
radiation. The transmission factors in turn determine the shape and magnitude of the 
power spectrum. Approximate expressions for 7a were recently presented in the literature 
[19]. In this study, we determine them numerically by solving the radial equation (2.12). 
Our numerical method is outlined in section 3.2. 
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3. Numerical methods 



In this section we briefly describe the numerical methods we employ for computing the 
angular eigenvalues and the transmission factors. 

3.1 Angular eigenvalues 

The allowed values of the angular momentum quantum numbers I, j and m are restricted 
[29] . The azimuthal number m can take any integer value, and I and j are positive or zero 
integers, such that 

l>j + \m\ and I - j - \m\ = 2k, where k G {0, Z+}. (3.1) 

The second equation in (3.1) implies that only even [odd] values of I are allowed when 
(j + |m|) is even [odd]. Note that the angular equation (2.11) has two regular singular 
points at 9 = and 9 = ir, just like the 4-D spin-weighted spheroidal equation (equation 
(2.11) with j = n = 0), plus an extra regular singular point at ^ = tt/2, which is a non- 
singular point in the 4-D case. The presence of this extra regular singular point means 
that half the regular solutions in the 4-D case will now, for n > 0, become irregular. This 
is the reason why the angular momentum quantum number I is forced to have the same 
parity as j + (as opposed to being allowed to be any non-negative integer in the 4-D 
case). Note also that the whole of the physical region in the bulk is covered by the regime 
9 G [0, 7r/2]. However, when restricting ourselves to the brane (by taking 9i = 0), the whole 
physical region is covered by extending the regime to ^ G [0, tt] (in the same way as when 
vertically slicing a 2-sphere into two halves we must extend the regime of the polar angle 
from 9 G [0, tt] to G [0, 27r] in order to cover the full circle). 

The numerical method we used for solving the angular equation in the bulk, equation 
(2.11), is essentially the same - with the obvious modifications - as the shooting method 
for the 4-D case as described in detail in [41], to which we refer the interested reader. We 
require power series expansions around the two integration end-points ^ = and 9 = 7r/2. 
Near the end-point 9 = 7r/2 we used the power series expansion given in equation (3.4) of 
[29]. Near the other end-point = (or = tt, by symmetry) we insert the expansion 

oo 

5a = (sin 0)1"! (cos 9y ^ (1 ± cos 9f (3.2) 

p=0 

into equation (2.11) and obtain the recursion relation: 

Oip 
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where 

Up = 2(p+ |m| + 1), 

Pp = [3{p - 1) + 2{n + 2j) + 4(|m| + l)]p- (aw)^ - + j{j + n + 1) 

+ \m\ (n + \m\ + 2j + 1), 
7p = 3{aujf -p{p + n + 2j + 2\m\ -l)+n + EA- j{j + n - 1) (3.4) 

— \m\ (n + \m\ + 2j — 1), 
Sp = — 3(au;)^, 
ep = (aw)^. 

It is easy to check that forn = j = this recursion relation coincides with the corresponding 
one in the 4-D case (see, e.g., equations (2.5)-(2.8) with spin s = in [29]). 

The eigenvahics E\ in the 4-D case become degenerate for large acu [41]. Fortunately, 
this degeneracy does not occur in the bulk case [29]. Therefore, the integration of the 
angular equation with our technique is more amenable in the bulk than in 4-D for large 

OLO. 

3.2 Transmission factors 

To compute the transmission factors 7a we solved the radial equation (2.12) numerically. 
We started with the ingoing solution close to the horizon and integrated out to large r^- 
Here, we matched the numerical solution onto the asymptotic form (2.22a) to determine 

and A'^^ and hence, via (2.34), T\. Typically, we integrated from = 1.001 to a 
matching point at = 100. 

The solution close to the horizon may be found using the method of Probenius. The 
radial equation (2.12) may be rewritten 

cPRa A{ri) dRA Biij) 

—nr = 1 — I — 2~ y^-^i 

drj-' ri dr] ry^ 
where 77 = r/r^ — 1. Here, A{r}) and B{rj) are functions with Maclaurin series expansions. 

The coefficients Aj and Bj are easily determined using a mathematical software package 
such as Maple. Near the horizon, the radial solution may be expressed as a power series in 
V, 

j=0 

The index s satisfies the indicial equation s{s — 1) — sAq — Bq = 0, which has solutions 
s = -iziuj / {2k) . The negative root is the correct choice for an ingoing solution. The 
coefficients Rj are determined from the recurrence relation 

[(s + j){s + j-l)-{s+ j)Aa - 5o] Rj = XI + i - k)Ak + Bk] Rj-k- (3.8) 

k=l 
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We used this method to compute expansion coefficients up to R^. 

The ingoing and outgoing solutions in the large-r^ Hmit may be found by a similar 
method. First we make the substitution 

^AW = ^W^^A(r), neZ+. (3.9) 

The positive (negative) sign in the exponent gives the outgoing (ingoing) solution. This 
substitution leads to a differential equation which can be written 

d^XA (±2iud ^ C{z)\ dXA D{z) 

where z = l/r^. Again, C{z) and D{z) have Maclaurin series expansions which can be easily 
found with a symbolic algebra package. The radial function Xf^^^z) also has a Maclaurin 
series expansion, X\{z) = 'Ylfj=Q^j^'' ^ with coefficients Xj determined by the recurrence 
relation 

i 

±2iud{j + l)X,+i = j{j - l)Xj - ^((j - k)Ck + Dk)Xj^k. (3.11) 

fc=0 

We computed coefficients Xj up to tenth order. Note that the ingoing solution is simply 
the complex conjugate of the outgoing solution. 



4. Results 



In this section we present numerical results for the emission of scalars by a rotating higher- 
dimensional black hole. Our objective is to compare the relative magnitudes of emission 
into the bulk and on to the brane. First, we begin with a review of the results for the 
Schwarzschild (non-rotating) phase [12]. Next, we examine the effect of black hole rota- 
tion upon the transmission factors and emission spectra. In its rotating phase, the black 
hole is shedding both mass and angular momentum; we determine the loss rates for both 
quantities. We conclude by presenting bulk-to-brane emission ratios for a range of values 
of and n. 

4.1 Non-rotating black holes (a^ = 0) 

A detailed study of scalar brane and bulk emission in the non-rotating phase was presented 
in [12]. Here we review some key results of that work. 

Figure 1 shows the power spectrum on the brane (left) and in the bulk (right) for 
n = 1, 2 ... 6. The total emission rate, or power, is equivalent to the area under the curves. 
The total emission depends strongly on the number of bulk dimensions, n. For all n, the 
total power emitted on the brane exceeds that emitted into the bulk (note the y-scales in 
figure 1). This is despite the fact that particles emitted into the bulk are, on average, more 
energetic than those emitted on the brane. 

Table 1 shows the total power (in units of 1/r^) for n = 1, 2 ... 6. The total power 
increases monotonically with n. The proportion of the total power that is emitted into the 
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Emission on tlie Brane Emission in tlie Buil^ 




Figure 1: Brane and hulk emission from a non-rotating higher- dimensional black hole. The left plot 
shows the power emitted on the brane, for various numbers of space-time dimensions {n = 1, 2 ... 6). 
The right plot shows the power emitted into the bulk. Note the factor of 3 difference in the scales 
on the y-axis. 





n = 1 n = 2 n = 3 n = A n = 5 n = 6 


Total Power, rldE/dt 
% in bulk 


0.00371 0.0134 0.0362 0.0823 0.170 0.339 
28.3% 19.9% 17.9% 19.6% 24.8% 34.0% 



Table 1: Proportion of scalar power emitted into the bulk by a non-rotating black hole. These 
quantities were calculated by numerically integrating the power spectra up to a frequency cutoff of 
ujrh — 10. The values are in good agreement with earlier studies [12]. 

bulk (i.e. the 'missing' energy fraction) decreases for n from 1 to 3, but increases again for 
n > 3. 

Figure 2 shows how the proportion of energy lost in the bulk depends on n. The total 
power emitted in the bulk appears to obey a power law relationship, rf^^ ~ Aa"', where 
A ~ 0.0004 and a ~ 2.54. The power emitted on the brane increases less rapidly with re. 
Hence the proportion of total energy lost in the bulk decreases up to re = 3, but increases 
for re > 3. 

The key question that we seek to address in the next section is: how does the proportion 
of energy lost into the bulk change when the black hole is rotating? In other words, is the 
standard claim [22] that black holes radiate mainly on the brane still correct in the spin- 
down phase of black hole evolution? 

4.2 Rotating black holes, > 

Here we consider a black hole rotating in the = 7r/2 plane on the brane. Introducing 
a single plane of rotation has many consequences: it modifies the angular eigenvalues, 
the transmission factors, and the form of the mode sum. Rotation breaks the azimuthal 
degeneracy, so that modes of different m contribute different amounts of power. 
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Total Power Emitted in Brane and Bulk Channels (a = 0) 
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Figure 2: Total power (in scalar particles) emitted by a non-rotating black hole. The blue line 
shows the power emitted in the bulk, and the green line shows the power emitted on to the brane, 
for various n (the number of extra dimensions). Note the logarithmic scale on the y-axis. 

Previous studies [15, 18] have shown that, overall, rotation enhances the emission on 
the brane. Higher angular momentum modes become increasingly important, and the black 
hole emits particles at higher energies. The superradiant effect increases the emission from 
the co-rotating modes, particularly for low n. For fast-rotating black holes, the m = I 
modes dominate the emission. This leads to an oscillatory power spectrum, and a rapid 
loss of black hole angular momentum. 

4.2.1 Transmission factors 

Let us begin by considering the transmission factors T\, defined by equation (2.34). An 
analytic expression for the transmission factor was derived in [19] under the assumption 
that both the energy uj^ and angular momentum parameters are much lower than unity. 
For the indicative case of the I = j = m = mode, the dominant behaviour, as a; ^ 0, was 
also computed and found to be given by the expression 



where A* = {n + 1) + {n — l)a^ and = 1 — Aa^/Al. In figure 3, we compare the 
approximation of equation (4.1) with our numerically-determined transmission coefficients. 
We find good agreement in the regime T\ <0.1. It is clear that the effect of rotation is to 
shift the I = j = m = transmission curve to lower energies. 





A, 2"(n + l)r2(Mi)(2-D,)' 
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Figure 3: Transmission factor of the I ^ j ~ Q mode for n — b, for a^i — 0.0 and Od — 1.2. This 
plot compares the approximation [19] (dotted Une) given in equation (4.1) with our exact numerical 
results (solid line). The right plot shows the same data, but with a logarithmic scale on the y-axis. 

To compute the total cross section, it is not sufficient to consider only the I = j = m = 
modes. We must compute a sum over all angular modes. The modes are labelled by 
three integers: /, the total angular momentum; m, the angular momentum in the plane of 
rotation; and j, the angular momentum in the bulk. These integers satisfy the condition 
(3.1). If the black hole is rotating then we must compute transmission factors for each valid 
combination of I, j and m. In practice, we will truncate the mode sum at some suitable I = 
^max- Unsurprisingly, the multiplicity of modes in the bulk is greater than the multiplicity 
of modes on the brane. For example, the number of non-degenerate modes N{L) in the 
range < / < L is A''branc(-^^) = {L + 1)^ on the brane, but is A'buik(-^) = L{L'^ — l)/6 
in the bulk. The computation time increases commensurately with the number of modes 
considered. 

Figure 4 shows the transmission coefficients for the first few modes, for the particular 
case of a slowly-rotating 6D hole (a^ = 0.4). A hierarchy of effects is clear. The modes are 
split firstly on I, secondly on m and finally on j. The azimuthal splitting on m increases 
in magnitude with rotation parameter a^. 

In the frequency regime w < rriQh, the superradiance effect causes the transmission 
coefficient to become negative. Nonetheless, the emitted power is positive, since the denom- 
inator of equation (2.32) is also negative in this regime. Superradiance is most significant 
for the m = I, j = modes (i.e. the maximally-corotating modes). Figure 5 shows how su- 
perradiance affects the m = I transmission factors of fast-rotating 5D (left) and 6D (right) 
black holes. It is clear that superradiance is most significant at low n. In particular, for 
n = 1, superradiance has a significant effect on the shape of the power spectra, as we see 
in the next section. 
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Figure 4: Transmission factors for the first few modes, I — 0, 1 & 2. The plot shows transmission 
factors for a 6D black hole rotating at = 0.4. The three sets of lines correspond to / = (left), 
I = 1 (middle) and I = 2 (right). For each I, there is a multiplet of (/ + !)(/ + 2)/2 modes. The 
modes in each multiplet are split by m and, to a lesser extent, by j. 
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Figure 5: Transmission factors for the maximally-rotating m = I modes of a fast-rotating (od = 
1.5^ black hole. The left plot shows transmission factors for a 5D (n = 1) black hole, and the right 
plot shows the same modes for a 6D (n = 2) black hole. The transmission factors are negative for 
uj < mQ.h due to the superradiance effect. 

4.2.2 Power spectra 

Let us now consider the effect of rotation on the shape and magnitude of the spectra of 
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scalar radiation emitted in both the bulk and the brane. As we saw in the previous section, 
the shape and magnitude of the emission spectrum strongly depends on the number of extra 
dimensions, n. Let us examine two cases in detail, n = 2 and n = 6. 



Emission on the Brane 




Emission in the Bull^ 




Figure 6: Brane and bulk emission from a 6D rotating black hole. The left and right plots show 
the power emitted on the brane and in the bulk, respectively. Black hole rotation increases the 
proportion of the total flux that is emitted on the brane. Note the order-of-magnitude difference in 
the scales on the y-axis. 

Figure 6 shows how rotation changes the emission spectrum of a 6D (n = 2) black 
hole. The left panel shows emission on the brane, and the right panel shows emission in 
the bulk. From the left panel, it can be seen that the total emission on the brane (the 
area under the curve) increases monotonically with the rotation parameter a^. In the 
bulk, the situation is more complicated. A slowly-rotating black hole actually emits less 
power into the bulk than a non-rotating hole. However, above ~ 0.4, the total emission 
increases with a^. The total powers for the curves shown in figure 6 are given in table 
2. The reduction in bulk emission is due to the fact that the Hawking temperature (2.30) 
decreases monotonically with ad- Thus, a rotating black hole is cooler than a non-rotating 
hole. On the brane, the reduction in Th is more than outweighed by the fact that rotation 
changes the transmission factors, and enhances emission in the higher angular modes. In 
the bulk, this is not necessarily the case, particularly at low ad- 





ad = 0.0 ad = 0.2 ad = 0.4 ad = 0.6 ad = 0.8 ad = 1.0 


Bulk Power 
Brane Power 
% in Bulk 


0.00267 0.00262 0.00256 0.00266 0.00308 0.00391 
0.01074 0.01148 0.01414 0.02003 0.03132 0.05142 
19.9% 18.6% 15.3% 11.7% 9.0% 7.1% 



Table 2: Proportion of scalar power emitted into the bulk by a 6D rotating black hole. These 
figures were calculated by numerically integrating the power spectra, up to a cutoff of ujrh ~ 10. 
The power (mass loss rate) is given in units of 1/r^. 



-16- 



In the non-rotating case, emission is primarily in the / = and I = 1 modes. However, 
rotation enhances the emission from higher angular modes. The oscillatory pattern in the 
emission spectra at higher ad is due to the increasing importance of the superradiant m = I, 
j = modes. This is seen clearly in figure 7, which shows the contribution of the m = I 
and m = I — 1 modes to the overall power spectrum. 



0.0014 



0.0012 



0.001 



T3 

ci^ 0.0008 



0.0006 



0.0004 



0.0002 



n = 2, a = ^.5r^ 
bulk 



4 5 



Figure 7: Emission from a 6D black hole at Od = 1.5. The plot shows that the majority of the 
emission comes from the maximally-corotating modes (m = /). The m — I — 1 modes are also 
plotted. 

Next, let us consider emission from a lOD (n = 6) rotating black hole. Figure 8 
compares the power spectra on the brane (left panel) and in the bulk (right panel), at 
Ud = 0.0 (red), ad = 0.5 (green) and = 1.0 (blue). It is immediately clear that rotation 
enhances emission on the brane significantly, but has a lesser effect on emission in the 
bulk. A large number of angular modes contribute to the total power for n = 6. Figure 
9 compares the relative magnitudes of the m = l, m = l — 1 ... ,m = I — 5 modes, at 
Ud = 1.5, for I = 0, 1, ... 26. Note that a sum over j was taken where multiple modes 
exist with the same / and m. With six bulk dimensions, the degeneracy factor Nj (2.35) 
appearing in the mode sum is potentially very large. This enhances the contribution from 
the higher j modes. The degeneracy factors increase in the sequence 1, n + 1, n(n + 3)/2, 
(n + 5)(n + l)n/6, ... for j = 0, 1, 2, 3, . . .; i.e. 1, 7, 18, 77, ... for n = 6. Modes close 
to m = have a much higher degeneracy than modes near m = I. For each m, the most 
significant contributions arise from the j = I — m modes. 

Figure 10 compares the emission on the brane (left) and in the bulk (right) at ad = 1.0, 
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Figure 8: Brane and bulk emission from a lOD rotating black hole. The left and right plots show 
the power emitted on the brane and in the bulk, respectively. Black hole rotation increases the 
proportion of the total power that is emitted on the brane. Note the order-of-magnitude difference 
in the scales on the y-axis. 

for a range of n = 1, 2 ... 6. Note that the scale on the y-axes differs by a factor of ~ 7. 
Comparing these spectra with the non-rotating spectra of figure 1, it is clear that rotation 
has reduced the proportion of the overall emission which enters the bulk. Table 3 lists the 
fraction of power lost in the bulk, for a range of scenarios. This data is plotted in figure 
11, which shows the % of scalar power in the bulk for = . . . 1 and n = 1, 2, ... 6. 
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ad = 0.5 


ad = 1.0 


n 
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28.3 
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0.00371 


20.9% of 


0.00449 
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0.01154 
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= 2 


19.9 


% of 


0.0134 


13.5% of 


0.01953 


7.1% of 


0.05533 


n 


= 3 


17.9 


% of 


0.0362 


11.8% of 


0.05497 


6.2% of 


0.1646 


n 


= 4 


19.6 


% of 


0.0823 


13.0% of 


0.1275 


6.8% of 


0.3808 


n 


= 5 


24.8 


% of 


0.170 


16.7% of 


0.2609 


9.1% of 


0.7709 


n 


= 6 


34.0 


% of 


0.339 


24.0% of 


0.5041 


14.7% of 


1.3832 



Table 3: Proportion of total power emitted into the bulk. This table gives the percentage of the 
total power (i.e. the combined bulk and brane mass loss rate, in units of l/r^) that is emitted 
into the bulk, for a range of dimensionalities n, and angular momenta Od- The total powers were 
calculated by numerically integrating the power spectra up to cutoffs of urh = 10, 12 and 16, for 
ad — 0.0, 0.5 and 1.0, respectively. 

Finally, let us briefly consider very fast rotation, a^^ > 1. As we have seen, a large 
number of modes are excited by a fast-rotating black hole. For instance, figure 9 shows 
that modes up to / ~ 30 contribute to the emission spectrum of a lOD hole at = 1.5, and 
that the power spectrum is significant at ur^ ~ 15 and beyond. Such spectra may give a 
misleading impression, since the black hole cannot emit a pair of particles which are more 
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Figure 9: Emission from a lOD black hole at ad = 1.5. The upper plot shows the total power. 
The lower plot shows the contribution from a 'forest' of modes, m = ^ to m = ^ — 5. It can be seen 
that the m = I modes are no longer dominant at all energies. Note that a sum over j has been 
taken where appropriate. 



energetic than the black hole mass itself. For black holes created at a collider, we would 
expect the black hole mass to be only a few multiples of M^f. Many other physical effects 
may modify the spectrum for uorh ^ 1. A concerted effort is currently underway to model 
the particle showers produced by black hole decay using Monte Carlo event generators [42] . 

Regardless of these physical considerations, it is still possible to compute an 'idealised' 
semi-classical emission spectra using equation (2.32). The calculation is easiest at low n, 
since many subdominant modes (m < / — Tn-cutoff) can be neglected. Figure 12 shows bulk 
emission from a 5D black hole at = 0, 1, 1.5, 2 and 3. Again, the oscillatory structure 
is due to the dominance of the superradiant {m = I) modes. In every case we have studied 
we find the power spectrum reaches an overall peak and tends to zero as wr/^ oo. 

4.2.3 Angular momentum spectra 

A rotating black hole loses both mass and angular momentum through the Hawking emis- 
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Figure 10: Emission by a rotating black hole. These plots compare the power emitted on the 
brane (left) and in the bulk (right) at — 1.0. Note the difference in the scales on the y-axis. 
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Figure 11: Total scalar emission for a range of scenarios. The upper plot shows the total power 
in brane and bulk scalars, for Od = 0, 0.5 and 1, and for a range of numbers of extra dimensions, 
n = 1, 2 ... 6. The lower plot shows the % of this power which is emitted into the bulk. 



sion process. The loss rate of angular momentum is given by equation (2.33). 

Figure 13 compares the loss rate of angular momentum on the brane (left panel) and in 
the bulk (right panel), for a range of rotations of a 6D (n = 2) black hole. Note the order- 
of-magnitude difference in their magnitudes. Hence, the majority of angular momentum is 



-20- 




0123456789 

Figure 12: Emission in the bulk from a fast-rotating 5D black hole. The plot shows how the power 
emitted into the bulk varies with black hole rotation rate, Od = a/rh = 0, 1, 1.5, 2 and 3, with one 
extra dimension n = 1. 

emitted on the brane. This is perhaps not surprising, since the plane of rotation hes on the 
brane. Table 4 hsts the total angular momentum emission rates for n = 2, ad = . . . 1.0. 
The proportion of angular momentum entering the bulk is found to decrease monotonically 
with ad- 
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Table 4: Proportion of (scalar) angular momentum emitted into the bulk by a 6D rotating black 
hole. These figures were calculated by numerically integrating the angular momentum spectra, up 
to a cutoff of Lorh = 10. The angular momentum loss rate is given in units of l/r/j. 

Figure 14 shows the angular momentum emitted by a rotating black hole embedded 
in a ten-dimensional bulk. Comparing with figure 13, we see that the overall emission rate 
has been enhanced by an order-of-magnitude by the increase in bulk dimensionality. We 
estimated the total angular momentum loss rate by integrating these spectra up to a cutoff 
of OT/j = 16. We found that, at ad = 0.5, some 9.0% of the total angular momentum loss 
rate of 0.28r^^ is emitted in the bulk. At = 1.0, we estimate that 6.0% of a total of 
1.64r^^ is emitted in the bulk. Note that the latter estimate is an underestimate due to 
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Figure 13: Brane and bulk angular momentum from a 6D rotating black hole. The left and right 
plots show the angular momentum emitted on the brane and in the bulk, for a range of angular 
momenta ad = a/rh- Note the order-of-magnitude difference in the scales on the y-axis. 



the integration cutoff (see figure 14). 
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Figure 14: Brane and bulk angular momentum from a lOD rotating black hole. The left and right 
plots show the angular momentum emitted on the brane and in the bulk, respectively. Note the 
order-of-magnitude difference in the scales on the y-axis. 



4.3 Consistency with other studies 

In this section we briefly compare our results with those published in the literature. We 
find our resuhs are consistent with the previous studies, where comparison is possible. In 
ah the papers mentioned below, the methods used to compute the greybody factors are 
different from those employed here, and therefore represent an independent check on our 
calculations. 

First, let us consider the emission of scalars in the Schwarzschild phase. We find 
excellent agreement with the results presented in [12], for both brane and bulk emission, 
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although it should be emphasized that the numerical algorithms used in [12] are different 
from those we have used in section 3. The bulk-to-brane proportions in table 1 agree well 
with those in table 6 in [12] (to ~ 1). We also find excellent agreement with the emission 
ratios presented in table 1 of the second paper in [13]. 

Next, let us consider scalar emission on the brane in the rotating phase; here, we are in 
agreement with [15], which again uses different numerical methodology to that in section 
3. The total powers on the brane presented in table 3 are in reasonable agreement with 
the totals in figure 12b in [15]. Our values are likely to be more accurate than those shown 
in figure 12b of [15], because we use a higher cutoff frequency {ujrh ~ 16 vs. uirh ~ 3). 

We have checked the credibility of our new results for scalar emission in the bulk in a 
variety of ways. For example, in section 4 we confirmed that the numerical greybody factors 
agree with the analytic approximation derived in [19] in the appropriate limit {ujrh < /). 
Comparison was also made with two numerical studies for the special cases n = I and 
n = 2 (i. e. 5D and 6D black holes). A study of bulk emission from a 6D rotating black 
hole lying on a tense brane was recently conducted in [21]. We confirmed with the authors 
that, when the brane tension is set to zero, our power spectra are in excellent agreement 
with their results (for the special cases n = 2, = 0, 0.6 and 1.2). On the other hand, 
we cannot reproduce the results of a recent study of a 5D hole [20]; whilst we find good 
agreement for brane emission, wc cannot verify their results for bulk emission (figure 5, 
[20]). Since this appears to be the only point of disagreement between our study and the 
existing literature, we remain confident that our numerical results are substantively correct. 

5. Discussion and conclusions 

The emission of Hawking radiation, in the form of elementary particles, by a higher- 
dimensional black hole might be the most distinctive signature of the creation of a black 
hole, and thus of the existence of extra dimensions in nature. Whereas the emission on 
the brane may in principle be directly observable and can be used as a source of valuable 
information about the extra dimensions, the emission in the bulk - the space transverse to 
our brane - will only be interpreted as a missing energy signal. Obviously, the more energy 
the black hole emits in the bulk, the less remains for emission on the brane. On the other 
hand, the amount of missing energy, due to the emission in the bulk and the emission of 
elusive particles such as neutrinos on the brane, is now considered as another distinctive 
signature of a black hole event. For the above reasons, the study of the emission in the 
bulk by a higher-dimensional black hole is of equal importance as the one on the brane. 

The bulk emission channel for a higher-dimensional black hole in its spherically- 
symmetric Schwarzschild phase was the first one to be studied in the literature. Two 
types of particles are usually considered as being allowed to propagate in the bulk: gravi- 
tons and scalars. A comprehensive study for the case of the latter type of particles [12] 
demonstrated that the emission in the bulk remains always subdominant to the one on the 
brane, although it may become important for a large number of extra spacelike dimensions. 
A number of studies in the case of gravitons [13] have finally concluded that, although the 
emission of gravitons is enhanced by the existence of extra dimensions more than any other 
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species of particles, the total emission of gravitational modes remains subdominant to the 
one corresponding to the Standard Model degrees of freedom on the brane. 

The study of the bulk emission during the spin-down phase of the life of the black 
hole has up to now been restricted to scalar fields. This is due to the increased complexity 
of the gravitational background around the rotating higher-dimensional black hole, and 
to the cumbersome formalism necessary to derive the equations of motion of higher-spin 
fields in the bulk. The existing scalar analyses themselves, although very illuminating, are 
either approximate [19] or refer to black holes with a specific dimensionality [20, 21]. For 
this reason, here we have attempted to present a comprehensive study of the bulk emission 
of scalar fields by a rotating black hole with arbitrary angular momentum and living in a 
spacetime with arbitrary number of extra dimensions. 

The quantization of scalar fields propagating in a higher-dimensional, rotating black- 
hole background closely follows the one for their 4-dimensional analogues, therefore, in 
section 2 we have presented only the main steps of this analysis. It is worth mentioning 
that well-known techniques [33, 34, 35, 36, 37, 38, 40] developed for 4-dimensional space- 
time are found to hold also in this case leading to similar expressions for the energy and 
angular-momentum emission rates. Our numerical methods, used to compute the angular 
eigenvalues and the transmission factors, were described in detail in section 3. Finally, 
oTir exact numerical results for the bulk scalar emission from a rotating black hole were 
presented in section 4. 

As a consistency check, our numerical analysis first reproduced the previously derived 
results for the emission of scalar fields by a non-rotating black hole in the bulk [12]. We 
then proceeded to study the rotating case, and demonstrated that the approximate anal- 
ysis presented in [19] was in very good agreement with the exact numerical one in the 
low-energy regime. Inevitably, deviations appear beyond this regime, necessitating a full 
numerical study. The computation of the value of the transmission factor at all energies for 
a rotating black hole demands the consideration of a very large number of angular modes, 
that increases as the energy increases too. Our final results for the transmission factor re- 
vealed a hierarchical splitting of the modes according to the angular momentum numbers 
(first on I, then on m, and finally on j) as well as the expected superradiance effect, that 
was most effective for the {m = I, j = 0) modes and for low values of n. Compared to 
the non-rotating case, the transmission factor for the most-dominant modes is found to be 
enhanced with the angular momentum of the black hole. 

The corresponding power spectra for scalar fields were then computed. For the purpose 
of comparing bulk and brane emission, we derived both types of spectra for the same set 
of values of a and n. We find that whereas the brane scalar emission is enhanced with 
the angular momentum of the black hole, in accordance with previous studies, the bulk 
emission is initially suppressed and starts increasing only after a certain value of a. This 
effect persists for all values of n and is caused by the fact that - contrary to the emission 
on the brane - the enhancement of the transmission factor in the bulk is not in a position 
to compensate for the decrease in the black-hole temperature with a. The complete picture 
emerges from the entries of table 3, which lists the proportion of the total power emitted 
by the black hole in the bulk. From these, we see that, for all examined values of a 
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and n, the bulk scalar emission is only a fraction of the brane emission. This fraction 
becomes important only for large values of n, as in the non-rotating case, and even then 
it is suppressed as a increases. We may thus safely conclude that, as far as the scalar 
channel is concerned, a rotating black hole emits a smaller fraction of its mass into the 
bulk than a non-rotating one. Our results therefore seem to add extra support to the 
argument presented in [22] that black holes radiate mainly on the brane. 

Finally, we considered the angular-momentum loss rate of the black hole. Perhaps not 
surprisingly since the black hole's sole angular momentum component lies on the brane, 
we found again that the black hole loses its angular momentum mainly by emission on 
the brane rather than in the bulk. Thus, for the indicative case of a 6-dimensional black 
hole we found that the proportion of the angular momentum emitted in the bulk can be 
described with single figures, as can be seen in table 4. 

The analysis of the brane-to-bulk energy balance is, however, not yet complete. The 
graviton emission channel for the rotating phase needs to be investigated before a final 
conclusion can be drawn. In the non-rotating case [13], it was found that the behaviour 
of the transmission factors for gravitons and scalar fields in the bulk shared a number of 
qualitative features, with the significant enhancement of the graviton power rate caused 
mainly by the larger number of degrees of freedom in the bulk. It remains to be seen 
whether, for a rotating black hole, the graviton transmission factors will exhibit the same 
characteristics as the scalar transmission factors found in this work. If so, it also remains to 
be seen whether the number of graviton states can tilt the balance towards bulk emission. 
For this, we would need the perturbed equation of motion for gravitons in a general rotating, 
higher-dimensional black-hole background which is unfortunately still missing from the 
literature. 
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